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Introduction 

In the paper [|]], Dunkl introduced the following difference-differential oper- 
ators acting on functions on a Euclidean space V, related to arbitrary finite 
groups G generated by orthogonal reflections in V: 

V e = ^+ £ (1) 

aeR+ [a,x) 

Here <9g denotes the partial derivative in direction £ £ V, R is the root system 
of the group G, i.e., the set of unit normals to the reflection hyperplanes, R + 
is its positive part with respect to some generic linear form on V, k a = k(a) 
is a C7-invariant function on R, s a is the reflection corresponding to the root 
a e R, and s a is the operator on the space of functions on V: 

Saf(x) = f(s a (x)). 

To be precise, Dunkl used slightly different operators, which are conjugated 
to (HD by the operator of multiplication by Y[(a,x) ka . 

The main property of the Dunkl operators is given by the following 
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Theorem 1 (Dunkl) The operators ([I]) commute with each other: 

[V 6 V,] = 0, (2) 

for all £, i] G V . 

The goal of this work is to describe certain generalizations of the Dunkl 
operators (H), preserving the property @. Some of these results were an- 
nounced in M. 

In Section 1 we consider generalizations of the form 

where A + is the set of unit normals to some set S of hyperplanes in V passing 
through the origin, A + is its positive part, and k a = k(a) is some function 
on A + . We show that the commutativity of the operators implies that S 
is the set of reflection hyperplanes of some Coxeter group G, A + = R + , and 
k is G- invariant. 

In the Section 2, we consider operators of the form 

V « = + ( a , Of x))s a , (4) 

where f a {z) are functions of one variable, not identically 0. 

The commutation relations (Q) are equivalent to a system of functional 
equations for the functions f a , a G R + . In the case when G is the Weyl 
group W of a simple Lie algebra, we show that, with the exception of A\, 
B-2, the only ^-invariant solutions, i.e., such that 

for all a G R + , are Dunkl's solutions, in accordance with 0. 

While the A\ case is trivial, the B2 case leads to the classical theory 
of Landen's and Jacobi's transformations of elliptic functions. We give the 
general VT-invariant solution in this case. This result has an interesting 
topological application in the theory of elliptic genera (see ||). 

In Section 3, we give a solution of the functional equations in terms of 
elliptic functions, for arbitrary reduced root systems. This solution gives 
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families of commuting differential-difference operators that we call elliptic 
Dunkl operators. 

We then show in Section 4, using techniques from [[J, ||, that these are 
essentially all solutions in the A n _i case. 

In Section 5, quantum elliptic Dunkl operators are introduced. These 
are pairwise commutative families of difference operators. They depend on 
a parameter /i, and are such that elliptic Dunkl operators appear in the 
first order term (semiclassical approximation) of their expansion in power of 
yU. These operators are related to the transfer matrices associated with the 



We conclude by discussing the possible applications of our results to the 
theory of integrable n-body systems. 

1. Operators of Dunkl type and Coxeter groups 

Let S be a finite set of hyperplanes in a Euclidean space V passing through 
the origin, A the set of unit normals to the hyperplanes in S (two for each 
hyperplane), A + the positive half of A with respect to some linear form on 
V, and k a some non-zero coefficients. 
Let be the operator 



Theorem 2 The operator and Vr, commute for arbitrary £ and r\ &V if 
and only if A + coincides with R + for some Coxeter group G and k(a) = k a 
is a G-invariant function. 

Proof: The commutator [Vf , can be rewritten in the form I+II+III (see 



i?-matrix of 0, §. 
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(5) 



U), where 



I 




II 



s a ) = 0, 



k a kp(a,/3)(:. 



1 



III 



■v 



(a,x)(s a ((3),x) 



S a S/3, 
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where (a,(3)^ v = (a, £)(/?, rj) — (a,r))(f3,£). One sees that to have a can- 
cellation of the terms in the sum corresponding to fixed rotation s a S/3, it is 
necessary that, together with a, (3, the vector 7 = s a ([3), also belongs to A. 
This means that A is the root system R for some Coxeter group. 

To prove that k is G-invariant, let us recall that if it is so, then [V^, V^] = 
according to Dunkl's theorem. Suppose that k is not invariant, i.e., there 
exist two roots a, f3 such that kp 7^ fc 7 , 7 = s a (/3). Subtracting from III the 
Dunkl identity with appropriate coefficients, gives a relation of the form 

where c a> p = 0, c 7jQ 7^ 0. But this is impossible, since the pole at (7,2;) = 
cannot be canceled. □ 

Thus we have shown that Coxeter groups arise naturally as a commuta- 
tivity condition for a natural generalization of Dunkl operators. 



2. Functional equations for the coefficients of generalized Dunkl 
operators 

Let us now consider for given Coxeter group G the following generalizations 
of Dunkl's operators: 

V« = ^+ ]T (a,0fa((a,x))s a . (6) 

aeR+ 

It is convenient to extend the definition of f a to all a G R, by setting 

/_«(*) = -fa(-z). (7) 

With this definition, is independent of the choice of the positive part R + 
of R. In fact, we can replace R+ by any subset of R consisting of one normal 
vector for each hyperplane, without changing V^. 

We may choose as before R to consist of unit vectors, but, since Weyl 
groups will be considered later, it is convenient to allow normal vectors to 
have arbitrary length, still preserving the condition that we have two nor- 
mal vectors a and —a for each hyperplane in S. At this point, this is no 
generalization, since the operators do not change if we replace a by a mul- 
tiple a' = ca and replace the corresponding function f a by the function 

fa>(z) = C _1 /a(c _1 2;). 
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A calculation analogous to the previous one leads to the following formula 
[V£,V„]= ( a ,P)t,vf<x(( a > x ))fp(( s M, x ))s a S/3, (8) 

a,/3£R+ 

where, as above, {a,P)^ v = (a, £)(/?, 77) — (a, t/)(/3, £). 
Theorem 3 The commutativity condition for the operators (|j) 

[V € ,V„] 

is equivalent to the system of functional equations 

Y, (a,f3)f a ((a,x))f p {(sM,x)) = 0, (9) 

a,/3£R+:s a sp=r 

for any given rotation r. 

Here (a, (3) denotes the oriented area of the parallelogram with sides a, (3, 
with respect to some orientation of the plane spanned by these two vectors, 
which is perpendicular to the rotation axis of r. Obviously, the equation @ 
is independent of the choice of orientation. 

Example: For the root system of type A 2 (see Fig. |l]), we have the following 
functional equation for the three functions associated with the three roots 
labeled in Fig. [I]. The functions associated to the other roots are then deter- 
mined by (0). 

f(x - y)g(x - z) + g(y - z)h(y - x) + h(z - x)f(z - y) = 0. (10) 

Let us now discuss G- invariant solutions of the functional equations fllPf ) : 
fg(a) = fa for all g G G, or, equivalently, 

= V g @g, for all g EG. 

We assume that none of the functions f a vanishes identically. We restrict 
ourselves to the case where G is the Weyl group W a semisimple Lie algebra. 
Obviously it is sufficient to consider the case of a simple Lie algebra. Ex- 
cluding the one dimensional case A\ where commutativity does not give any 
restriction on f a , the first case we consider is the case A n _i, where G = W 
is the symmetric group S n . 
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Figure 1: The root system of type A2 



Proposition 4 For root systems of type A n _i, n > 3, the only S n -invariant 
solution of (JD meromorphic in a neighborhood of the origin is 

fa(z) = — , 

which corresponds to the usual Dunkl operators. 

Proof: Consider the functional equation (g). For given r, it involves only 
roots in the 2-plane orthogonal to the rotation axis of r. This plane is spanned 
by any two roots a, (3 with s a sp = r. 

In the A n _i case, the roots in any such plane build a root system of type 
A\ x A\ or A 2 . In the former case roots are orthogonal and (|9|) is identically 
satisfied for all /. It is therefore sufficient to consider the case n — 1 = 2. Let 
us rewrite ([H]) in the form 

f(u)g(u + v) + g(v)h(-u) + h(-u - v)f(-v) = 0. (11) 

S^-invariant solutions correspond to / = g — h, where / is an odd function, 
since f a (z) = f-a( z ) — ~fa(—z). This leads to the following equation for /: 

f{u)f(u + v) - f{v)f(u) + f{u + v)f(v) = 0. 

If / vanishes at a point u, then f(v)f(u + v) = for all v, and / vanishes 
identically. If / does not vanish anywhere, we get, after substitution (ft = 1/ /, 

(f)(u + v)= (j)(u) + 4>(v), 
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which has the only solution <j>(z) = cz. □ 

Let us consider some more examples of low rank. It will be shown below 
that these examples essentially cover the whole theory. 

Example: In the case of G 2 (see Fig. 0), one has two A 2 systems. Hence the 
the invariant solutions have the form 

f\z) = — > 9{z) = — , 

z z> 

for arbitrary A and B. 




Figure 2: The root system of type G 2 



Example: In the B 2 case (see Fig. |3|) one has the following functional equation 
for the symmetric solution: 

f{x){g{x + y)+g(x-y)) + f{y){g{x + y)-g{x-y)) = 0, (12) 

where / and g are odd functions. In this case we have more complicated 
solutions, such as 

f(z) = cot(z) g{z) 



sin (2) 

The general solution will be given below. 
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9 f 9 




Figure 3: The root system of type B 2 



Example: In the B 3 case, R = {±e», ±e, ± e^}, where {ex, e 2 , e 3 } is an 
orthonormal basis of V. Let us find the most general invariant solution. The 
roots a = t\ + e 2 , and (3 = e 2 + form a subsystem isomorphic to A 2 . Thus 
/« (^) = fp( z ) — C/z. Now consider the B 2 system generated by 7 = e\ and 
a. So we have to solve the functional equation ( Jl2| ) with g = C/z: 

f,(x)(-^- + — ) + / 7 (y)(-^ — ) = 0, 

or, equivalently (if C ^ 0), 

z/rOs) ~ yfiiv) = °> 

which implies / 7 (^) = D/z. It follows that the most general invariant solution 
is Dunkl's solution. 



Theorem 5 For all root systems of simple Lie algebras except Ax, B 2 , all 
W -invariant solutions of the functional equation (|9|) have the form 



fa(z) = — , 



where a ^ k a is a W -invariant function on R. 



Proof: For all root systems except Ax, B 2} there exists a subsystem isomor- 
phic to A 2 (see ||10|| ) All other two-dimensional subsystems are isomorphic 
either to A 2 or B 2 (the exceptional case of G 2 was considered above). Con- 
tinuing as in the B3 case, we complete the proof. □ 
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In the rest of this section we consider the invariant B2 case. In other 
words, we want to find the general solution of the equation (|12| ) 



f(x)(g(x + y)+ g(x - y)) + f(y)(g(x + y)- g(x - y)) = 0, (13) 

which are meromorphic in the vicinity of the origin. 

There are the following obvious solutions: f(z) = 0, g an arbitrary func- 
tion, and g(z) = 0, / an arbitrary function. We call these solutions trivial. 

Theorem 6 The general non-trivial solution of the functional equation (|13|) 
has the form 

or, more explicitly 

g(z)= , A 7a , /(*) = B ~. , (14) 

sn(az,k) sn(eaz,k) 

where k = (1 — k)/(l + k), e = — i(l + k), and A, B and a are arbitrary 
constants. 

Here sn(z,k) is the classical Jacobi elliptic function (see, e.g., JllJ). In the 
degenerate cases we have the following solutions: 

az 
cot — , 
2 ' 

B 



/(*) 


A 

smaz ' 




M 


= A cot az 






A 

z 


9{z) = - 



sm az 



with arbitrary constants A, B, a. In particular, we see that all non-trivial 
solutions are odd functions, and thus lead to invariant Dunkl operators for 
root systems of type B 2 . 

We proceed to prove the Theorem. 
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Lemma 7 If f(z), g(z) satisfy the equation (|13|) ; then the same is true for 
f(z), g{z), where 

(i) /O) = g(z) = ng(az), (15) 

for arbitrary constants X, \i, a, or 

(ii) f{z)=g{z), ~g{z)=f{z/2). (16) 

Proof: The first symmetry is evident. To prove the second one, it is sufficient 
to change variables 

u = x + y, v = x — y. 
Then equation fll3|) takes the form 



9(u ,(/(^) + /(^)) +9M (/(^)-/(^))^o. 

□ 

We will use these symmetries again and again, in particular in the proof 
of the following 

Lemma 8 There are no non-trivial solutions of ([13D, for which f(z) or g(z) 
is regular at the origin. 

Proof: Because of the symmetry (16), it is enough to consider only the 



case when f(z) is regular at the origin. Putting y = in (0), we get 
2f(x)g(x) = 0, which means that the solution is trivial. □ 



Lemma 9 Non-trivial solutions of ([13]) are odd functions: 
f(-z) = -f(z), g(-z) = -g(z). 



Proof: Rewrite the equation (13) in the form 



<t>{y) {g{x + y)+ #0 - v)) + <P(x) {g{x + y)- g{x - y)) = 0, 

where 4>(z) = l/f(z) is regular at vanishes at z = 0. By putting x = in 
this relation, we obtain 

<!>(yMv) + 9(-v)) = o 
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which implies for a non-trivial solution 

g(-y) = -g(y)- 

The fact that / is odd follows then from the symmetry (|TB]) . □ 
Let us introduce A = 0'(O), where (f>(z) = 1/ f(z), as before. 

Lemma 10 If f(z), g(z) is a non-trivial solution of (|13|) , then 

(i) A = 0'(O)/O, 

(ii) g = "*/(*), (17) 

H A/(* + y) = P{x) _ p{y) • ( 18 ) 

Proof: Rewrite equation (|H|) in the form 

/ , \ f(y) - IM ( \ / 10 n 

g(g + y) = /(y) + /(x) g(a; - y) - (19) 
Taking the logarithm of both sides and applying the operator d x + d y , gives 
g'(x + y) f(x)f'(y)-f(y)f(x) 



(20) 
(21) 



g(x + y) f 2 (y)-f 2 (x) 

<j)'(x)(j)(y) - <p'(y)<p(x) 
cj) 2 (x) - <p 2 (y) 

Putting y = in (]2~lD implies 

In particular, if A = 0, g is constant and therefore regular, which is impossi- 
ble. Now the formula (|18[) follows from (II7D and fl20T) □ 



Using the symmetry (|T5D, ([[§), we may set without loss of generality 
A = 1, so that 

f(z) = - + 0(z), g(z) = - + 0(z), (*-0). 
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The function / satisfies the functional equation (addition theorem) 

fix I v) ~ mm - mf(x) (22) 
Rewrite it in the following form 

and expand the right hand side near y — 0, using <j>{z) = z + az 3 + 0(z 5 ): 

fix + y) = f(x) + f\x)y + (3a f(x) + f{x))y 2 + 0(y 3 ). 

By comparing this with the Taylor expansion, one has 

f"(x) = 2f(x) + 6af(x), (24) 

implying, after multiplication by /' and integration, that 

[f'f = f 4 + 6af 2 + b, (25) 

for some constant b. The function <ft = 1/ f is thus a regular odd solution of 
the equation 

(0') 2 = 1 + 6a0 2 + 60 4 , (26) 
and therefore coincides with the Jacobi elliptic function 

<Kx) = Sn{€X e ' k \ (27) 

where (1 + k 2 )e 2 = —6a, k 2 e A = b. Recall that sn is the solution of the 
equation (s') 2 = (1 — s 2 )(l — k 2 s 2 ), with initial condition s(0) = 0. It satisfy 
the addition formula discovered by A. Cayley (see [PH ) 

s 2 (x) - s 2 (y) 



s(x + y) 



s(x)s'(y) - s(y)s'(x) 



This implies the relation ( p2]) for f(x) = e/sn(ex, k). Note that by ( JIB]) , also 
g has the same form as /, in general with different values of e and k. 
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Lemma 11 If f(x) = e/sn(ex,k) and g'(x)/g(x) = —f(x) then f(x), g(x) 
are a non-trivial solution of the functional equation 



Proof: We have 

g'(x + y) 



-2f(x + y) 



g(x + y) 

n /'(*)/(v) - /'WW 



So 



f 2 ( x ) - P(y) 
f(y) - fix) m + fix) 

m - m m + m ' 

f(y) - f(x) 



(d x + d y ) \ogg(x + y) = (d x + d y ) log 



f(y) + f( x Y 

This means that 

/ . \ f(y) - f( x ) i { \ 
fw + f( x ) 

for some function if). Rewriting this equation as 

g(x + y) = if){x-y), = 1//, 

(p{x) + (j)(y) 

and putting y = 0, gives ip{x) — g(x), and therefore 

/ \ f(y) ~ f( x ) , \ 
f(y) + f{ x ) 

which is equivalent to fllS]). □ 

To finish the proof of Theorem ^, we have to prove that the relation 
between k, k, e in 

1 ^ 

9( x ) = —, — TT' f( x ) 



sn(x, k) ' ' sn(ex, k) 

where g'/g = —f is the one given in Theorem |6]. 

The function 1/ f = —g/g' = sn(ex, k)/e satisfies the equation 

[(g/g')'} 2 = [i - t\g/g r m - ~k 2 e 2 (g/ 9 n (28) 
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or, equivalently, 



(g') 2 - gg") 2 = ((g'r - eV) ((</) 2 -^V 

Substituting (g'f = (g 2 - l)(g 2 - k 2 ), g" = 2g 3 - (k 2 + l)g, into (||), gives 
the following two possibilities 

(i) 2k = -(k 2 + 1) - e 2 , -2k = -{k 2 + 1) - k 2 e 2 

(ii) -2k = -{k 2 + 1) - e 2 , 2k = -{k 2 + 1) - k 2 e 2 

In the first case, we get k = (1 — k)/(l + k), e — i(k + 1), as required. The 
second case leads to an equivalent answer. Theorem |6] is proved. 

Remark. The transformation of elliptic functions g — > f is of second order 
and therefore can be reduced to the well-known Landen transformation (see 
[p"T|j). One can check that it is the composition of the "imaginary Jacobi 
transformation", which is the unimodular transformation with the action on 
the homology of the curve described by the matrix J 



1 
-1 



J = 

and Landen's transformation with matrix L 



2 
1 



Let us call it LJ-transformation. Thus the functional equation (|T3| ) describes 
a pair of elliptic functions related by the LJ-transformation. This fact has 
found recently an interesting topological application in the theory of elliptic 
genera M. 

Note that / and g, as elliptic functions of second order "live" on different 
elliptic curves, one of which is a double cover of the other. In particular, the 
solutions found here are not special cases of the elliptic solutions considered 
in the next Section. 
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3. Elliptic Dunkl operators 

In this section we consider only the case where R is the root system of a 
semisimple Lie algebra with Weyl group G = W. Let us consider the elliptic 
curve with modular parameter t, Im(r) > 0, and the family of functions 

ax{z) = O&M-X) ' AeC \ Z + ^ Z > ( 29 ) 
given in terms of Jacobi's theta function 



0x{z) = - E 



e 



27ri(2+i)(n+i)+7riT(n+i) 2 



n=— oo 



The functions <J\ have the following defining properties: 

(i) *x(z+l) = (T X {Z). 

(ii) a x {z + r) =e 2mX a x {z). 

(iii) a\ is meromorphic, its poles are on the lattice Z + rZ, and cr\( z ) = 
1/2 + 0(1) as z -> 0. 

More properties of this functions are given in the Appendix. 

Theorem 12 For any generic A G Vc = V^rC, and W -invariant function 
k a , the functions 

f a (z) = k a a {a v iX) (z), 
where a y = 2a/(a,a), satisfy the functional equations (|9]). 

Proof: Fix the rotation r. All roots involved in the left hand side of the 
functional equation (0) 

I( x ) = H («,/5)o"(Q v ,A)((a,x))a(^v iA) ((s a (/3),x)) (30) 

a,/3£R+:s a S0=r 

lie on the same two-dimensional plane. Consider J(x) as a meromorphic 
function of x G Vc, and let P v = {p G V|(p, a) G Z Va G R}. Then if 
p G P v , /(s + p) = I(x), and as a; — > x + pr, the term labeled by (a, (5) in 
the sum (^0]) gets multiplied by 

e 27ri(( a v ,A)(a,p)+(/3 v ,A)( SQ (/3),p))_ 
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Since s a (X) = A — (a v , A) a, we see that the multiplier can be rewritten as 

c ) 

and is therefore the same for all terms in the sum. It follows that I(x) has 
the quasi-periodicity property 

I(x + q + pr) = e 2m{x - rW ' p) I(x), q + pr G P v + rP v . (31) 

Let us now consider the poles of the function I. Poles may appear when x is 
on the hyperplanes (a,x) = 0, a G R + , or their translates by P v + rP v . If 
x approaches the hyperplane (a, x) = 0, the singular terms in the sum (j30|) 
are the terms indexed by a,/3 and 7,5 where s 7 (5) = ±a. As all roots are 
on a plane this implies that 7 = =F/?. The sign is + (thus s 1 (S) = —a) since 
7 > 0. In particular only two terms are singular in (pCp. The coefficient of 
the (simple) pole is (cf. (iii) above) 

(a,(3)a^v jX )((s a (P),x)) - ((3,5)a {f3 v :X) ((f3,x)). 

This expression vanishes on the hyperplane (a,x) = because s a (x) = x 
there, and 

{(3,5) = -(s p (P) lS p(6)) = -((3,a) = (a,f3). 

It follows that the singularity at (a, x) = (and thus on all affine hyperplanes 
(a,x) = n + rriT, n, m G Z by (EH])) is removable. We conclude that / has no 



singularity on Vc, and has the quasi-periodicity property (^T|). It therefore 
vanishes, by Fourier series theory. □ 

Corollary 13 The operators 

V\ = d(.+ k a (a,£)o-( a v tX) ((a,x))s a , 

form a commutative family: 

[v 5 A ,v, A ] = o. 

Let us call these operators elliptic Dunkl operators. They are not VT-invariant 
but VT-equivariant 

^V^- 1 = V:g. (32) 
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Example: In the A n _\ case, we identify functions on V = {R n | = 0} 
with functions / on R n such that f{x\ + a, . . . ,x n + a) is independent of 
a G R. Let ei, . . . , e n be the standard basis of R n and put 

1 n 

e« = ei y~] ej G V. 

71 — 

« j = l 

Then Dunkl operators are linear combinations of the commuting operators 

V, A = — + fc a*-^ ( x i ~ ( 33 ) 
OXi jW* 

where s^- is the operator that interchanges the ith and jth variable. 

Remark. It is interesting to note that the usual Dunkl operator reminds 
Moser's L-matrix for the Calogero problem whereas the elliptic Dunkl 



operator (|3^) reminds Krichever's generalization [DJ, see also |T^]. The gen- 



eral solution of the corresponding functional equation (different from ours) 



was found by Bruschi and Calogero [15]. More general functional equations 



motivated by these problems, as well as some topological problems, were 
introduced and solved by one of the authors 0, @. 

4. General solution of the functional equation in the A n ^i case 

In this section we show that elliptic Dunkl operators are essentially the only 
solutions of the functional equation (|9]), in the A n _i case, n — 1 > 2. 



As before, we start with A 2 . In this case, the functional equation is (|TTD 



f(u)g(u + v)+ g(v)h(-u) + h(-u - v)f(-v) = 0. (34) 

We must find the most general solutions of this functional equation, where /, 
g, and h are assumed to be meromorphic functions defined in a neighborhood 
of the origin. 

First of all, if one of the three functions vanishes identically, then the 
functional equation says that the product of the other two vanishes, and we 
get a solution where two functions vanish and the third is arbitrary. We call 
these solution trivial, and consider from now on only solutions where /, g 
and h are not identically zero. 
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Lemma 14 Iff, g, h satisfy Q3|), then (i) f(z) = g(z), g(z) = h(z), h(z) = 
f{z), (ii) f(z) = af{bz)e az , g(z) = ag(bz)e l3z , h(z) = af(bz)e yz ! for arbitrary 
constants a, . . . , 7 such that a + (3 + 7 = 0, 
also satisfy fl34]) . 



Proof: Replacing by (— u — v,u) in (|3^) implies (i). Property (ii) is 

easy to check. □ 

Proposition 15 The only non-trivial solutions of ( p4[ ) holomorphic around 
the origin are 



f(u) = ae au , g{u) = be^, h{u) = ce 
where a, b, c ^ 0, ab + be + ac = 0, a + (3 + 7 = 0. 

Proof: It is easy to check that these are solutions. We prove uniqueness. 
Suppose /, g, h are a non-trivial solution of (|34"1), defined and holomorphic 
in a neighborhood of the origin. Taking v = in (|34|), we get 



/(t*M«) + (/(0) + </(0))/i(-«) = 0. (35) 

We have /(0) + ^(0) 7^ since does not vanish identically. Similarly, by 
Lemma [14|, 

g{u)h(u) + (g(0)+h(0))f(-u) = 0, (36) 

h(u)f(u) + (h(0)+f(0))g(-u) = 0. (37) 

Introduce the functions F(u) = f(u)f(—u), G(u) = g(u)g(—u), H{u) = 
h{u)h(-u), and S{u) = f{u)g{u)h{u). Multiplying (|§, (0) by h{u), 
f{u) and g'('u), respectively, we obtain 

S{u) = XF(u) = fiG(u) = uH(u), (38) 

for some non-zero constants A, fi, v. In particular S is an even function, and 
thus 

S(u) 2 = S{u)S{-u) = F{u)G{u)H{u). 

Hence F(u) 3 = constF(u), and F is a constant. Similarly, G and H are 
constant functions. Thus we have 



f{u)f(-u) = a 2 , g(u)g(-u) = b 2 , h(u)h(-u) = c 



2 
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for some constants a, 6, c 7^ 0. We can therefore write 

f{u) = ae^ u \ g{u) = be^ {u) , h{u) = ce^ u \ (39) 

for some odd functions 0, ip, 77, holomorphic in a neighborhood of 0. Since 
S(u) = f(u)g(u)h(u) is constant, we have + ip + r\ = 0. Inserting (^) in 
the functional equation gives 



ahe 4>{u)^iu+v) + hce i>{v)- v {u) + ace -„( u+ „)-0(„) = Q 



which, after elimination of rj = — — ip can be recast in the more convenient 
form 

ab + bc^(u,v) + ac$(w, v)' 1 = 0, (40) 

with $(w, v ) = exp(0(-u+t>)— 0(w)— 0(t>)), and ^(w, t>) = exp(ip(u+v)— ip{u) — 
i/j(v)). In particular, setting u = v = 0, we obtain the condition ab+bc+ac = 
0. Moreover, we have $(— w, — f ) = $(«, and it, — = ^(u, u) -1 , 
since and -0 are odd functions. Replacing (w,t>) by (— w, — 1>) in ( f4"0"D yields 
the equation 

ab + bc^(u, v)~ l + ac$(w, = 0. (41) 

Elimination of \1/ from (^), fl4"T| ) gives a non trivial quadratic equation with 
constant coefficients for $. Thus $ is constant, implying that \1/ is constant 
as well, we conclude that 0(w + f) = 4>{u) +0(t> ), and ip(u+v ) = ip(u)+ip(v), 
which leaves us with the solution 0(w) = au, ip(u) = j3u. □ 

From now on, we consider the case when f(u) has a pole of order p > 
at the origin. It is easy to see that h and g also have a pole of the same order 
p at u = 0. 

Inserting the expansion 

/(«) = £ + -, (42) 
9(u) = £ + (43) 
%) = 4 + ---> (44) 



into (|34]) , and multiplying the relation by u p (u + v) p v' p , yields a series in u 
and w which starts as 

acu p + abv p + (-l) p bc(u + v) p . 
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We see that p = 1 is the only case in which cancellation is possible. In this 

case, 

ac — be = 0, ab — be = 0, 

which implies a = b = c. Without loss of generality, we consider the case 
a = b = c = 1. 

Suppose that /, g, h are a solution of (|34f) with a simple pole at the origin 
with unit residue: 

f(u) = - + /o + ---, (45) 
g(u) = ~ + <7o + -'-, (46) 
h(u) = - + ho + - • • . (47) 

The left hand side of (0) has a Laurent expansion at v — 0: 

+ (~ + <k,)M-«) + (M-^) - + /o) + o(v). 

The constant term gives the equation 

f(u)g(u) + (f + go)h(-u) + h'{-u) = 0. (48) 

Similarly, by Lemma [14] (i), 

g(u)h(u) + (g + ho)f(-u) + f'(-u) = 0, (49) 
h(u)f(u) + (h + f )g(-u)+g'(-u) = 0. (50) 

Introduce F(u) = f(u)f(—u), G{u) = g(u)g(—u), H[u) = h(u)h(—u), and 
S(u) = f(u)g(u)h(u), as above. Note that F, G and H are even functions. 
Then F'(u) = f'(u)f(-u) - f'(-u)f(u), and © implies that F'(u) = 
S(u) — S(—u). More generally, (^8| -p0|) imply 

F\u) = G\u) = H\u) = S{u) - S(-u). 

We can then write 

F(u) = a- P(u), G(u) = b-P(u), H(u) = c- P(u), (51) 
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where P(u) = 1/u 2 + 0(u 2 ) has no constant term, and a, b, c are integration 
constants. The even function P obeys 



P'{u) = -S(u) + S(-u). 



(52) 



Let us compute the derivative S' = (f'/f + g' / g + h'/h)S using (pj H50|) 

S'(u) = aS{u)-G{u)H{u)-H{u)F{u)-F{u)G{u) 
= aS(u) -3P(u) 2 + /3P(u) +7, 

for some constants a, (3, 7. Let us consider the even and odd part of this 
equation separately: 



±(S(u)+S(-u)) 
d 



du 



(S(u) - S(-u)) 



a(S(u) - S(-u)) = -aP'{u) (53) 
a(S(u) + S(-u)) - QP(u) 2 + 2(3P(u) + 2 7 . (54) 



Integrating (|53|) gives 

S(u) + S(-u) = -aP(u) + 5, 
for some 5. By subtracting (|52|) from this equation, we obtain 



S(u) 



'-aP(u) 



P\u)). 



(55) 



(56) 



Inserting (|55D and (^) into (p4]) yields finally the differential equation 

P"(u) = QP(u) 2 + (a 2 - 2j3)P{u) - a5 - 2j. 

Let us multiply this equation by P'{u) and integrate. We get 

P'iu) 2 = 4P(u) 3 + CiP(u) 2 + c 2 P(u) + c 3 . 

In fact Ci = 0, since, by construction, the Laurent expansion of P has no 
constant term. This equation is well-known to have a unique meromorphic 
solution with double pole at the origin. It is the Weierstrass function p 
for some elliptic curve determined by the coefficients C2, C3, or one of its 
degenerations ^ 2 (1/ sin(a;w) 2 — 1/3), 1/w 2 (see [O]). After rescaling of the 
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variables as in Lemma O if necessary, we may assume that the periods of 



p are 1 and r, and that cu = n. By writing the constants a, b, c in (pT|) 
as -P(A), -P(/i), P{v) respectively (this is always possible since P defines a 
surjective map from the elliptic curve C/Z + rZ, or, in the degenerate cases, 
(C/Z)U{ioo}, CU{ioo}, onto the Riemann sphere), we see that the equation 
f(u)f(—u) = P(A) — P(u) has the solution (see the Appendix) 

f(u) = a x (u), 

and, in the degenerate cases, 

„, . ii sin(7r(w — A)) , , . , 

t(u) = -. — r -. — = 7r(cot(7ru) — cot(7rA)), 

v ' sm(vra) sm(-TrA) v v ; v " 

f(u) = i-I. 

u A 

The general solution in a neighborhood of the origin is f(u) exp <f)(u) for some 
odd function <p(u) regular at the origin. Similar formulas hold for g and h: 
the functions 



f(u) = a\(u), g{u) = a^u), h(u) = a v (u), 

and their degenerations, are a solution of the functional equation, provided 
that A + /i + v = (modulo the lattice). We now show that these are the 
only (up to transformations of Lemma |14l (ii)) solutions such that (|51|) holds 
with a = -P(A), b = P(fi), c = -P(^). Suppose /, g, h is another solution. 
Then 

/(«) = e*Wf(u), ~g{u) = e+Mg(u), h{u) = e^h(u), (57) 

for some odd functions 0, if), rj. Since the product S = fgh is expressed in 
terms of P (see flSED), the sum + ip + V vanishes identically. Inserting (|5"7j) 
in (^9|), (|50D, we deduce immediately that <p'(u) = ip'{u) = for all u. Thus 
4> and ip are linear functions. 

Combining this result with Prop. p~|, we obtain the following theorem. 

Theorem 16 The following list exhausts all non-trivial solutions of the func- 
tional equation 

f{u)g{u + v ) + g(v)h(—u) + h(— u — v)f(—v) = 0. 
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Elliptic solutions: 

f(u) = aa\(bu)e au , g(u) = aa ll (bu)e? u ', h(u) = aa v {bu)e lu \ 

where the function o~\ is defined in ( P^fl . 
Trigonometric solutions: 

aMHu-X)) 'y-f A (58) 

sin (em J sin(- oAj sm{bu) sm(— o//J 



a sin(6(n — i/)) 
sin(cw) sin(— eV) 



3 7« 



Rational solutions: 



f(u) = (-- i)e", = (" - V", %) = (- - V" (59) 

u A u jj u v 

Regular solutions: 

f(u) = -\e au , g(u) = --e^, h(u) = -~e^. (60) 

A H h> 

The parameters in these solutions are complex numbers a, b, a, j3, 7, A ; fi, 
v satisfying a + /3 + 7 = 0, A + /i + z/ = 0. 

In the trigonometric and rational cases, the limiting cases where X, fi, v 
take the value ±ioo are permitted. 

Remark. Notice that all these solutions are limiting cases (degenerations) of 
the elliptic solutions. 

This theorem covers the case A 2 , but gives immediately the answer in the 
case A n _i, n — l>2. In this case, we identify, as above, functions on V with 
translation invariant functions on R". 

Corollary 17 Let n > 3. The n operators 

Vj di -\- ^ ] fij(%i -Ej^Siji i 1, . . . , 71, 

with fij ^ 0, fij(u) = —fji(—u) are pairwise commutative iff fij(x) = 
kaxi-x ■(bx)e( ai ~ a ^ x where a\ is the function defined in (129]) , or one of its 
degenerations (see above). 
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5. Quantum Dunkl operators 

Fix complex parameters r, \x and k such that Im(r) > 0, \i ^ Z + rZ and 
k ^ 0. Consider the operator R(\), depending on the "spectral parameter" 
A G C, and acting on the space of, say, meromorphic functions of two complex 
variables X\ and x 2 : 

R(X)f(xi, X 2 ) = — 7Tt{0" m (xi2 + %)f(Xi + £, a7 2 - £) - 0a(x12 + ^)/(x 2 , X X )}. 

We use the notation xyi to denote the difference x\ — x 2 - The operator i? 
obeys the quantum Yang-Baxter equation 

i?(A 12 )( 12 ) J R(A 13 ) (13) i?(A 2 3)( 23 ) = J R(A 2 3) (23) J R(A 13 ) (13) J R(A 12 )( 12 ). 

The two sides of this equation are operators acting on functions of three 
variables, and the notation R(X)^ indicates the operator -R(A) acting on a 
function of several variables, by viewing it as a function of the ith and jth 
variable. This solution of the Yang-Baxter equation is essentially the one 
introduced in || as three-parameter generalization of the two-parameter so- 
lution of Shibukawa and Ueno [|7| . For positive integer values of k it admits a 
restriction to a finite dimensional subspace coinciding with Belavin's solution 

With the normalization used here we have "unitarity" 

i?(A 12 )( 12 )i?(A 21 ) (21) = Id, 

and R(0)f(xi, x 2 ) = f(x 2 , X\). Moreover -R(A) tends to the identity as /i goes 
to zero. 

Solutions of the quantum Yang-Baxter equation with these properties can 
be used to construct commuting operators Tj(Ai, . . . , A„) (related to transfer 
matrices of integrable models of statistical mechanics): 

Tj(Ai, . . . , A„) = R {i ' i+1) ■ ■ ■ R^R^ ■ ■ ■ R^- 1 ). (61) 

These operators act on functions of n variables, and, as before, R™' = _R(Aj — 
is the operator i?(Aj — A^) acting on a function of n variables by viewing 
it as a function of the ith and jth variable. 

Our result is that elliptic Dunkl operator defined in ( |33D can be ob- 
tained as semiclassical limit of the quantum operators, if k is integer. 
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Theorem 18 Let Tj(Ai, . . . , A n ) be the operators ( |61| ) acting on functions f 
on R n ; such that f(x\ + a, . . . , x n + a) = f(x\, . . . , x n ) for all a 6 R. For all 
integer k, we have 

Ti{\ x ,..., A n ) = Id + ^g-jy^n/n + 0(a* 2 ), 

where g m is the function 

9m(x) = n^i(^j) m exp(m^x^ / 1 (Ai i )/ei(Ai j )), 

viewed as multiplication operator, and the parameter k of R is given by k = 
(—l) k nk. 

Proof: Let p(x) = 8[(x) / 9 1 (x) . We first compute the expansion of R to first 
order. 

R(X) = Id + nr( A) + 0(/i 2 ), 
where r, the "classical r-matrix", is the differential-difference operator 

r W = -(^i - 9 2 ) + cr\(x 12 )si2 + p(A) - p(x 12 ). 

K 



This implies that 



T^Ax, . . . , A n ) = Id + fi £ r^\X i:i ) + 0(/x 2 ). 



Since Ti^dif = 0, we have 



and therefore 



X] (di -dj)f = ndif, 



Ti(Ai, . . . , X n ) = Id + -// ( di + ~~ <7 '\ (•''<./ + -(P(^ij) ~ P( x ij)) 1 • 



k \ rrz . n rrz . n 



The claim follows then from the relations 



dig m {x) = mJ2 (P(^ij) ~ P( x ij)) > 



Sij9m ( 1) 9m^ijj 

kl 



and the elementary fact that m = (— 1) k obeys (— l) m m = k. □ 
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6. Quantum n-body problems 

Trigonometric and rational Dunkl operators in the A n _i are used in the 
theory of integrable quantum n-body problems of the Calogero-Sutherland 
type, see |J, [J3J. Let us consider the trigonometric Dunkl operators (in a 
slightly different normalization) 

Vj = di + ^2 k(coth(xi - Xj) - 1)%. 

Then the operators Lj = ResX)"=i(Vi)- J , j — 1, . . . ,n form a set of pairwise 
commuting S^-invariant differential operators, and generate the algebra of 
S n -invariant differential operators commuting with the Schrodinger operator 

n i 
i=l #j Smh { X i - X j) 

Here Res M is the differential operator whose restriction to ^-invariant func- 
tions coincides with the differential-difference operator M. 

In the elliptic case, we are lead to consider the commuting operators 
Mj(X) = E?=i(^i) j , where V, A is the elliptic Dunkl operator fl33|). These 
operators are not ^-invariant. Instead, they obey wMj(X)w~ 1 = Mj(w(X)). 
However, let us consider the singular limit when A tends to the symmetric 
point 0, for j = 2. We use the shorthand notation 00 00 ^ 00 j • 

n 

M 2 (A) = £(V 4 A ) 2 
i=i 

n 

i=l ij^j ij^j 

In this calculation, we used the functional equation of Prop. |l^ (iii) for a\ 
and the fact that cr_ A (— x) = —a\(x). Since 

(r\{x)ax{-x) = p(X) ~ p(x), 

and limA-*o v'\{ x ) — ~p( x ) ~ 2t?i, (see the Appendix) we see that 

M 2 (A) - 5>(A«) 

i±3 
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has a limit as A — > 0: 

n 

L 2 = lim(M 2 (A) - ^2 p(\j) - ®i ~ + 1) P( Xi ~ x i) + const - 

v£j i=l v£j 

This is (minus) the Schrodinger operator of the so-called elliptic Calogero- 



Moser integrable n-body problem (see, e.g., pl|). 

It is reasonable to conjecture that the higher S'n-invariant differential 
operators commuting with L 2 can be obtained as A — > limits of suitable 
equivariant polynomials in with A dependent coefficients. 

For integer k there is a conjecture of one of the authors (see |L7|) that 
there are additional integrals of motion, such that the whole ring of quantum 
integral is supercomplete. We hope that elliptic Dunkl operators will help to 
prove it. 

Another interesting problem is to understand the analogue of Opdam's 



shift operator in the elliptic case (cf. [fig]). The one-dimensional case shows 



that it can not be a pure differential operator, because the genera of the 
spectral curves for Lame operators depend on the integer parameter k. 

In the B 2 case our results imply the quantum integrability of the system 
with Hamiltonian 



H 



where F = f — f 2 , G = g 



A + F(x) + F(y) + G(x + y) + G(x - y), 

g 2 , and /, g, are given by the formula (|I 



The commuting operator K has the form Res^fV^), where Vj are the cor- 
responding generalized Dunkl operators (0). The quantum integrability of 



this systems was independently established in [19 



It would be interesting to understand the relation between this construc- 
tion and the recent construction of elliptic Dunkl operators of [20] , which are 
formal infinite linear combinations of affine Weyl group reflections. 
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Appendix: Properties of the function cr\(z) 



We give some properties of the function a\(z) (j29|). This function can be 
expressed in terms of the Weierstrass functions (see, e.g., with periods 
1, r: 



a{z)a(—X) 2 



Proposition 19 (i) <t\{z) = —a-\(—z), 

(ii) cxaO) = -0i(A), 

(iii) a x {z)a fl (w) - cr x+fl (w)a x (z - w) - a^w - z)a x+lx (z) = 0, 

(iv) ax(z)a-x(z) = p(z) - p(X), 

(v) \im x ^ (T\(z) = -p(z) - 2?7i. 



Proof: (i) and (ii) are obvious. The functional equation (iii) follows from the 
fact the left hand side is 1-periodic in, say, z, and is multiplied by exp(27riA) 
as z — > z + r. Moreover, the singularity as z — ► (and thus as z — > m + nr, 
n, m G Z, by quasiperiodicity) is removable since (T\(z) = 1/z + 0(1). Thus 
we have an entire 1-periodic function of z with non-trivial multipliers as 
z — > z + r. Such a function must vanish. 

(iv) The left hand side is an even elliptic function with double pole with unit 
coefficient at the origin, and no other pole in the period parallelogram. Thus 
it is of the form p(z) + 4>(\), for some function if). Then use cr\(\) = to 
determine if). 

(v) Recall that a'(z) = ((z)a(z), ('(z) = -p{z) and a{z) = z + 0(z 5 ). We 
have 

d a{z-X) a(z-X) {({z - A) - ({z)) 

dza(z)a(—X) a(z)a(—\) 
= C'(z)+0(X). 

Thus c' x (z) = —p(z) + 0(A) + 2r]i\(T\(z), but \a\(z) — > — 1 as A tends to 0. 
□ 

Remark. Note that the functional equation (iii) is a special case (A2) of 



Theorem 12 
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